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Abstract 

This is a primer on extended Gaussian Images. Extended Gaussian Images are 
useful for representing the shapes of surfaces. They can be computed easily from: 

1. Needle maps obtained using photometric stereo, or 

2. Depth maps generated by ranging devices or stereo. 

Importantly, they can also be determined simply from geometric models of the 
objects. Extended Gaussian images can be of use in at least two of the tasks facing 
a machine vision system: 

1. Recognition, and 

2. Determining the attitude in space of an object. 

Here, the extended Gaussian image is defined and some of its properties discussed. 
An elaboration for non-convex objects is presented and several examples are shown. 
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1. Introduction 


In order to recognize an object and to determine its attitude in space, it is 
necessary to have a way of representing the shape of its surface. Giving the distance 
to the surface along parallel rays on a regularly spaced grid provides one way of 
doing this. This simple representation is called a depth map. A range finder produces 
surface descriptions in this form. Unfortunately, depth maps do not transform in a 
simple way when the object rotates (For one thing, interpolation must be used to 
get a new depth map on a regularly spaced grid). Alternatively, surface orientation 
might be given for points on the surface on some regular sampling grid. This 
grid may conveniently correspond to the picture cells in an image. Such a simple 
representation is called a needle map (Figure 1) [13]. Photometric stereo produces 
surface descriptions in this form [7, 12, 15, 28, 32, 33]. This representation too is 
not directly helpful when it comes to comparing surfaces of objects that may be 
rotated relative to one another (Both depth maps and needle maps depend on the 
position of the object as well as its attitude). 

The extended Gaussian image, on the other hand, does make it easy to deal 
with the varying attitude of an object in space [2, 3, 9, 13, 16, 29]. For one thing, 
it is insensitive to the position of the object. Some information appears to be 
discarded in the formation of the extended Gaussian image, however. Curiously, in 
the case of convex objects, the representation is nevertheless unique. That is, no 
two convex objects have the same extended Gaussian image. 

This representation of the shape of the surface of an object has proven most 
useful in work on automatic bin picking [4, 17, 18]. We start our discussion 
with objects having planar faces. Later we consider smoothly curved objects. 
Methods for computing discrete approximations of extended Gaussian images, 
called orientation histograms, are presented too. Orientation histograms can be 
computed from experimental data or -mathematical descriptions of the objects. 
Sections marked with an asterisk may be omitted on first reading or if your interest 
in the mathematical details is limited. 

2. Discrete Case: Convex Polyhedra 

Minkowski showed in 1897 that a convex polyhedron is fully specified (up to 
translation) by the area and orientation of its faces [21, 23, 26]. We can represent * 
area and orientation of the faces conveniently by point masses on a sphere. Imagine 
moving the unit surface normal of each face so that its tail is at the center of a unit 
sphere. The head of the unit normal then lies on the surface of the unit sphere. This 
sphere is called the Gaussian sphere and each point on it corresponds to a particular 
surface orientation. The extended Gaussian image of the polyhedron is obtained 
by placing a mass at each point equal to the surface area of the corresponding face 
(Figure 2). 

It seems at first, as if some information is lost in this mapping, since the 
position of the surface normals is discarded. Viewed another way, no note is 
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Figure 1. A needle map shows unit surface normals at points on the surface on 
a regular grid. Normals which point towards the viewer will be seen as dots, while 
tilted surface patches give rise to normals which are shown as lines pointing in the 
direction of steepest descent. 

made of the shape of the faces or their adjacency relationships. It can nevertheless 
be shown that (up to translation) the extended Gaussian image uniquely defines 
a convex polyhedron [29]. An iterative algorithm has recently been invented for 
recovering a convex polyhedron from its extended Gaussian image [22]. 

2.1. Properties of the Extended Gaussian Image 

The extended Gaussian image is not affected by translation of the object. 
Rotation of the object induces an equal rotation of the extended Gaussian image, 
since the unit surface normals rotate with the object. 

Mass distributions which lie entirely within one hemisphere, that is, are zero 
in the complementary hemisphere, do not correspond to closed objects. As we shall 
see, the center of mass of an extended Gaussian image has to lie at the origin. This 
is clearly not possible if a whole hemisphere is empty. Also, a mass distribution 
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Figure 2. The extended Gaussian image of a polyhedron can be thought off 
as a collection of point masses on the Gaussian sphere. Each mass is proportional 
fx • to the area of the corresponding face. Point masses on the visible hemisphere are 

shown as solid marks, the others as open marks. The center of mass (shown as the 
symbol 0) must be at the center of the sphere if the polyhedron is a closed object. 


which is non-zero only on a great circle of the sphere corresponds to the the limit 
of a sequence of cylindrical objects of increasing length and decreasing diameter 
(Figure 3). We will exclude such pathological cases and confine our attention to 
closed, bounded objects [23, 29]. 

Some properties of the extended Gaussian image are important: Firstly, the 
total mass of the extended Gaussian image is obviously just equal to the total 
surface area of the polyhedron. If the polyhedron is closed, it will have the same 
projected area when viewed from any pair of opposite directions. This allows us to 
compute the location of the center of mass of the extended Gaussian image. 

Imagine viewing a convex polyhedron from a great distance. Let the direction 
from the object towards the viewer be given by the unit vector v. A face, with unit 
normal swill be visible only if h ■ v > 0. Suppose that the surface area of this 
face is Oj. Due to foreshortening it will appear only as large as would a face of area 

(§t -y)Oi, 

normal to v (Figure 4). The total apparent area of the visible surface is 

^(v) = £ (s» • v) Oi, 

{* 1 Si-V>0> 
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when viewed from the direction v. The total apparent area of the visible surface 
when viewed from the opposite direction is 

A(—v) = - Z (s* • v) Oi, 

{i|Si-V<0} 

This should be the same, that is, A(v) = A{ —v). Consequently, 

Z ( § » Z °i 

all t Lall t 

where the sum now is over all faces of the object. This holds true for all view 
vectors, v, so we must have 

£ k Oi = o. 

an i 

That is, the center of mass of the extended Gaussian image is at the origin. 

An equivalent representation, called a spike model, is a collection of vectors 
each of which is parallel to one of the surface normals and of length equal to the 
area of the corresponding face. The result regarding the center of mass is equivalent 
to the statement that these vectors must form a closed chain when placed end to 
end (Figure 5). 
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Figure 4. A surface element appears smaller because of foreshortening. The 
apparant area is the true area times the cosine of the angle between the surface 
normal and the vector pointing towards the viewer. 




Figure 5. Vectors parallel to the normals of the faces of a polyhedron, and 
of length equal to the areas of the corresponding faces, form a closed chain when 
placed end to end. 


2.2. Reconstruction of a Tetrahedron (*) 

Faces that share a common edge are said to be adjacent. The masses on the 
Gaussian sphere corresponding to two adjacent faces are not always each others 
closest neighbors. Recovering a polyhedron from its extended Gaussian image is 
not easy in the general case [22], because it is hard to determine which faces are 
adjacent. Finding the actual offsets of each of the faces from the center of mass of 
the polyhedron is not as hard. 
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Figure 6. A tetrahedron with vertices A, B, C, and D. We are to find the 
distances of the faces from the center of mass, given the areas and surface normals 
of the faces. 


The structure of a tetrahedron, however, is very simple: Every face is adjacent 
to the other three. The shape of the tetrahedron is completely determined by the 
surface normals of the four faces, only the size of the tetrahedron remaining to be 
determined. In other words: There is only one degree of freedom left. Another way 
to look at it is to note that the four faces must have areas that place the center 
of mass of the extended Gaussian image at the origin, as we have just seen. This 
condition places three constraints on the four parameters. 

A A 

Let the given unit surface normals be a, b, c, and d, and the areas of the 
corresponding faces, A, B, C, and D (Figure 6). We have to determine the distances, 
o, b, c, and d, of these faces from the center of mass of the tetrahedron. From these 
distances we can, if desired, compute the positions of the vertices A, B, C, and D, 
simply by intersecting three of the planes at a time. The notation here is that the 
face opposite vertex A has area A and unit surface normal a, and so on. 

The perpendicular distance of the center of area of a triangle from one of 
the sides is one third the perpendicular distance of the vertex opposite that side. 
Similarly, in a tetrahedron, the distance from the center of mass to a particular face 
is equal to one quarter of the distance of the vertex opposite that face. We start by 
finding a formula for the distance of the face with area D, say, from the opposite 
vertex d. The desired distance, d, will be just a quarter of the result obtained in 
this fashion. The remaining three distances, a, b, and c, can then be computed 
using formulae obtained by cyclical permutation of the variables. 
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The position of the reconstructed tetrahedron is abritrary, since the extended 
Gaussian image is insensitive to translation. To make the result unique, we might 
place the center of mass at the origin. To reduce the size of the expressions to be 
manipulated here, however, it is convenient to move the tetrahedron so that one 
vertex, D, say, is at the origin. The distances of the faces from the center of mass 
are obviously not affected by this. 


Suppose for now that we know the locations of the vertices A, B, and C relative 
to D. We can then compute the directions of the six edges of the tetrahedron, 
by taking all of the distinct pairwise differences of the four vertex positions. Four 
surface normals can then be found by taking cross-products of these edge-direction 
vectors. We actually need only four of the edge vectors forming a closed circuit to 
do this. The results can then be normalized to obtain unit surface normals, 

B X C * C X A . A X B 
3 fB X C|’ jCXAf C_ |A X B|’ 

and, 


(A — C) X (B — A) _ AxB-fBxC-f CxA 
|(A - C) X (B - A)| |AxB + BxC + CxAf 

Now the perpendicular distance of the plane with area D from the origin can be 
found by taking the dot-product of any of the three vertices, A, B, and C with the 
unit normal d. Thus, 


4d = <iA = dB = dC = 


[ABC] 

|AXB + BXC + CXA|* 


The area of the facet opposite the origin is also easy to compute, 


D = |(A — C) X (B — A)| = ^|A X B -f B X C + C X A|. 

Our task is to express the offset d in terms of the area D and the given unit surface 
normals. The two formulae above do not allow us to do that directly, because we 
do not know what the value of [ABC] is. This quantity, by the way, is six times 
the volume, V, of the tetrahedron, or, 



V=i(4d)D=T[ABC], 

We proceed by considering the four distinct triple products of the four unit surface 
normals. First of all, 


fa b cl =_[ABCp_ 

1 1 |A X B[|B X C||C X A| ’ 

since [(x X y) (y X z)(z X x)] == [xyz] 2 . Then, by similar reasoning, 


abd] = 


[ABC] 2 

B X C||C X A||A XB-fBxC-fCxAj’ 
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since [xy (x -f-y 4" z )] — [xyzj. Formulae for [bed] and [cad] can be found by 
cyclical permutation of the variables. 


Multiplying the three formulae found this way together we get, 

[ABC] 6 


A A 


abd] [bcd][cad] 


(A X B) 2 (B X C) 2 (C X A) 2 1A XB-f BxC-f CX A| 3 ’ 


and so 


[abd][bcd][£ad] |ABC] 2 j 

[a be] 2 |AXB + BXC + CXA| lJl h 

So that finally, 

a j \J(2D)[k b d] [b c cl] [c ad] 

4d —_ - ~ * ~ ■. 

— [a b c] 

The other distances, a, b, and c, can be computed using similar formulae obtained 
by cyclical permutation of the variables. 


3. Continuous Case: Smoothly Curved Objects 

The ideas presented in the previous chapter can be extended to apply to 
smoothly curved surfaces. 

3.1. Gaussian Image 

One can associate a point on the Gaussian sphere with a given point on a 
surface by finding the point on the sphere which has the same surface normal 
(Figure T) [6, 11, 23]. Thus it is possible to map information associated with points 
on the surface onto points on the Gaussian sphere. There are no two points with 
the same surface normal in the case of a convex object which has positive Gaussian 
curvature everywhere. The mapping from the object to the Gaussian sphere in this 
case is invertible: Corresponding to each point on the Gaussian sphere there is a 
unique point on the surface (If the convex surface has patches with zero Gaussian 
curvature, curves or even areas on it may correspond to a single point on the 
Gaussian sphere). 

One useful property of the Gaussian image is that it rotates with the object. 
Consider two parallel surface normals, one on the object and the other on the 
Gaussian sphere. The two normals will remain parallel if the object and the Gaussian 
sphere are rotated in the same fashion. A rotation of the object thus corresponds 
' to an equal rotation of the Gaussian sphere. 

3.2. Gaussian Curvature 

Consider a small patch 60 on the object. Each point in this patch corresponds 
to a particular point on the Gaussian sphere. The patch 60 on the object maps 
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Figure 7. The Gaussian image of an object is obtained by associating with 
each point on its surface the point on the Gaussian sphere which has the same 
r' surface orientation. The mapping is invertible if the object has positive Gaussian 

curvature everywhere. 



into a patch, 6S say, on the Gaussian sphere (Figure 8). If the surface is strongly 
curved, the normals of points in the patch will point into a wide fan of directions. 
The corresponding points on the Gaussian sphere will be spread out. Conversely, if 
the surface is planar, the surface normals are parallel and map into a single point. 


These considerations suggest a suitable definition of curvature. The Gaussian 
curvature is defined to be equal to the limit of the ratio of the two areas as they 
tend to zero. That is, 


K = 


I ™ 77; 
<50— »o <50 


dS 

dO' 


From this differential relationship we can obtain two useful integrals. Consider first 
integrating K over a finite patch 0 on the object: 


U Kdo =U ds=s ' 

where S is the area of the corresponding patch on the Gaussian sphere. The 
expression on the left is called the integral curvature. This relationship allows one 
to deal with surfaces which have discontinuities in surface normal. 

Now consider instead integrating l/K over a patch S on the Gaussian sphere: 
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IJl/KdS = IJd° = 0 , 

where 0 is the area of the corresponding patch on the object. This relationship 
suggests the use of the inverse of the Gaussian curvature in the definition of the 
extended Gaussian image of a smoothly curved object, as we shall see. It also 
shows, by the way, that the integral of 1 /K over the Gaussian sphere equals the 
total area of the object. 

3.3. Alternate Definition of Gaussian Curvature 

Consider a plane which includes the surface normal at some point on a smooth 
surface. The surface cuts this plane along a curve called a normal section (Figure 
9) [6, 11, 20, 24, 26]. Let the curvature of the normal section be denoted by 
kn- Consider the one-parameter family of planes containing the surface normal. 
Suppose that 0 is the angle between a particular plane and a given reference plane. 
Then /cjv varies with 0 in a periodic fashion. In fact, if we measure 0 from the plane 
that gives maximum curvature, then it can be shown that 

kn( 6) — Ki cos 2 6 -J- K 2 sin 2 6, 

where is the maximum and K 2 is the minimum curvature. These two values 
of k are called the principal curvatures. The corresponding planes are called 
the principal planes. The two principal planes are orthogonal, provided that the 
principal curvatures are distinct (Figure 9). 
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Figure 9. Normal sections of the surface are made with planes which include 
the surface normal. The planes corresponding to the largest and smallest value of 
curvature are referred to as the principal planes. The Gaussian curvature is equal 
to the product of the largest and smallest values of curvature. 


It turns out that 

K = 

is equal to the Gaussian curvature introduced earlier. This is clearly zero for a 
plane. It is equal to 1 /R 2 for a spherical surface of radius R, since the curvature 
of any normal section is 1/R. 

* A ruled surface is one which can be generated by sweeping a line through 

space. A hyperboloid provides one example of such a surface. Developable surfaces 
are special cases of ruled surfaces [6, 11, 24]. Cylindrical and conical surfaces are 
examples of developable surfaces (Figure 10). On a developable surface at least one 
of the two principal curvatures is zero at all points. Consequently the Gaussian 
curvature is zero everywhere too. 

3.4. The Extended Gaussian Image 

We can define a mapping which associates the inverse of the Gaussian curvature 
at a point on the surface with the corresponding point on the Gaussian sphere. Let 
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Figure 10. A conical surface is an example of a developable surface. On it 
the Gaussian curvature is everywhere zero, because (at least) one of the principal 
curvatures is zero. 


u, and v be parameters used to identify points on the original surface. Similarly, let 
£ and r) be parameters used to identify points on the Gaussian sphere (These could 
be longitude and latitude, for example). Then we define the extended Gaussian 
image as 

where (£, rj) is the point on the Gaussian sphere with the same normal as the 
point (u,v) on the original surface. It can be shown that this mapping is unique 
(up to translation) for convex objects. That is, there is only one convex object 
corresponding to a particular extended Gaussian image [1, 26, 29]. The proof is 
unfortunately non-constructive and no direct method for recovering the object is 
known. 



3.5. Properties of the Extended Gaussian Image 

The center of mass of the extended Gaussian image of a smoothly curved 
object is at the origin. We show this in a way similar to that used earlier for 
extended Gaussian images of polyhedral objects. Consider viewing a convex object 
from a great distance. Let the direction from the object towards the viewer be 
given by the unit vector v. A surface patch with unit normal s will be visible only 
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if s • v > 0. Suppose its surface area is <50 (Figure 4). Due to foreshortening it will 
appear only as large as would a patch of area 

(s • v) <50, 

normal to v. Let H(v) be the unit hemisphere for which s-v > 0. Then the apparent 
area of the visible surface is 

= /h,v,/ G ( S 

when viewed from the direction v. The apparent area of the visible surface when 
viewed from the opposite direction is 

A (~v) = j H( _^ I G ( § ) ( s • —v) dS. 

This should be the same, that is, A(v) = A{ —v). Consequently, 

f s / G(i) (5 • v) dS = (/ s / C(8) 8 is) ■ v = 0, 

where the integral now is over the whole sphere S. This holds true for all view 
vectors, v, so we must have 

JjG(i)idS = o. 

That is, the center of mass of the extended Gaussian image is at the origin (This, 
by the way, is not a very helpful constraint in practice, since one usually only sees 
one side of the object). 

Another property of the extended Gaussian image is also easily demonstrated. 
The total mass of the extended Gaussian image equals the total surface area of the 
object. If one wishes to deal with objects of the same shape but differing size one 
may normalize the extended Gaussian image by dividing by the total mass. 

We can think of the extended Gaussian image in terms of mass density on 
the Gaussian sphere. It is possible then to deal in a consistent way with places 
on the surface where the Gaussian curvature is zero, using the integral of 1/K 
shown earlier. A planar region, for example, corresponds to a point mass. This in 
turn corresponds to an impulse function on the Gaussian sphere with magnitude 
proportional to the area of the planar region. 

It is tempting, by the way, to imagine that one can find the attitude of an 
object by lining up the principal axes of inertia of the observed extended Gaussian 
image and the one computed from the geometric model [29]. This would be rather 
straightforward, requiring only the calculation of the eigenvectors of a three by 
three inertia matrix. In practice, one typically has information only about the 
visible hemisphere and thus cannot compute the required first and second moments 
over the whole sphere. 

3.6. Objects that are not Convex 

Three things happen when the surface is non-convex: 




13 



Horn 


EGI 


1. The Gaussian curvature for some points will be negative. 

2. More than one point on the object will contribute to a given point on the 
Gaussian sphere. 

3. Parts of the object may be obscured by other parts. 

We chose to extend the definition of the extended Gaussian image in this case to 
be the sum of the absolute values of the inverses of the Gaussian curvature at all 
points having the same surface orientation, 

1 


G(e,i) = £ 




This definition is motivated by the method used to compute the extended Gaussian 
image in the discrete case, as we will see later. 

The above extension makes sense if there are a finite, or at most a countable, 
number of points on the surface with the same orientation. At times, however, all 
points on a curve or even an area on the surface have parallel surface normals. In 
this case we may use, 


C(A) = /„/*(»-•)*>. 


where h is a unit vector on the Gaussian sphere, while s is a unit vector on the 
surface of the object. The integration is over the whole surface of the object O and 
6 is the unit impulse function defined on a sphere. 

We can be more specific, if we let r(u, w) be a vector giving the point on the 
surface corresponding to the parameters u and v, then 

G(£, r)) cos rj — j J 6(£ — 0(u,v), rj — <j>(u, u))|r„ X r v | dudv, 

where 0{u, v ) and <f>(u, u) are the latitude and longitude of the point on the Gaussian 
sphere which has the same orientation as the surface does at the point (u, v). A 
planar region of area A will thus contribute an impulse of weight A to the extended 
Gaussian image, while a cylindrical region will give rise to an impulse wall along a 
great circle at right angles to the axis of the cylinder. The integral of the impulse 
wall will be equal to the area of the cylindrical region. 

Usually we think of the extended Gaussian image as a fixed entity associated 
with an object. In the case of non-convex objects we might want to alter the definition 
to include only those parts of the surface visible from a particular direction. This 
would make the (modified) Gaussian image dependent on the view-point. We avoid 
this potential complication here. 

3.7. Examples of Extended Gaussian Images (*) 

The extended Gaussian image of a sphere of radius R is 

C(f, n) = R\ 
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Figure 11. Ellipsoid with contours obtained by cutting the surface with three 
orthogonal planes passing through pairs of points where the Gaussian curvature 
has stationary values. 


as discussed already. 

Perhaps slightly more interesting is the case of an ellipsoid with semi-axes a, 
b and c lined up with the coordinate axes (Figure 11). An equation for its surface 
can be written 

(;/+GO’+(;)■-• 

More useful for our purposes here is a parametric form 

x =a cos 9 cos (/>, 
y —b sin 9 cos <j>, 
z—c sin<j>. 

A normal at the point 

r = (a cos 9 cos <f>, b sin 9 cos <j>, c sin <f>) T 
on the surface is given by 


n = (be cos 9 cos (f>, ca sin 9 cos <p, ab sin <j>) T , 
as will be shown later. The Gaussian curvature turns out to be equal to 


K 


abc 

2 

abc' 

(be cos 9 cos <j>y -f- (ca sin 9 cos <j >) 2 -f- (ab sin 0) 2 


n? 


where n 2 


n ■ n. 
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Figure 12. Latitude and longitude can be used to identify points on the Gaussian 
sphere. Each point on the Gaussian sphere corresponds to a unique surface 
orientation. 


If we let £ be the longitude and rj be the latitude on the Gaussian sphere, then 
a unit normal at the point (£, r?) on the sphere is given by (Figure 12) 

n = (cos £ cos rj, sin £ cos rj, sin rj) T . 

Now n = nn. Identifying terms in the two expressions for surface normals at 
corresponding points on the ellipsoid and the Gaussian sphere we get, 

be cos 9 cos <f> — n cos £ cos rj, 
ca sin 6 cos </> = n sin £ cos rj, 
ab sin <f> = rising, 

so that 

n 2 [(acos £ cos?/) 2 -j- (b cos £ sin rj) 2 -f- (csin rj) 2 ] = ( abc ) 2 , 
and finally, substituting for n 2 in the equation for K we get 

abc 2 

(a cos £ cos r;) 2 -j- (6 sin £ cos r}) 2 -(- (c sin r]) 2 

The extended Gaussian image, in this case, varies smoothly and has the stationary 
values 





at the points where r is equal to (±1,0, 0 ) T , (0, ±1, O) 71 , and (C, 0, ±l) r , respectively. 
These results can be easily checked by sectioning the ellipsoid using the xy, yz 
and 2 z-planes. The Gaussian curvature, in this case, equals the product of the 
curvatures of the resulting ellipses. One then uses the fact that the maximum and 
minimum curvatures of an ellipse with semi-axes a and b are a/b 2 and b/a 2 . 

Later we will derive the extended Gaussian image of a torus, an object that is 
not convex. 


4. Discrete Approximation: Needle Maps 

Consider the surface broken up into small patches of equal area. Let there 
be p patches per unit area. Erect a surface normal on each patch. Consider the 
polyhedral object formed by the intersection of the tangent planes perpendicular 
to these surface normals. It approximates the original surface. The smaller the 
patches, the better the approximation. 

The extended Gaussian image of the original (smoothly curved) convex object 
is approximated by impulses corresponding to the small patches. The magnitude of 
each impulse is about 1/p, corresponding to the area of the patch it rests on (Figure 
13). Strongly curved areas will distribute their impulses over a large region on the 
Gaussian sphere, while areas which are nearly planar will have them concentrated 
in a small region. In fact, the number of impulses per unit area on the Gaussian 
sphere approaches p times the absolute value of the Gaussian curvature as we make 
p larger and larger. This can be shown using the integral of 1/K given earlier. 

The tesselation of the surface can be based on an arbitrary division into 
triangular patches as long as the magnitude of each impulse on the Gaussian 
sphere is made proportional to the area of the corresponding patch on the surface. 
Alternatively, one can divide the surface up according to the division of the image 
into picture cells. In this case one has to take into account that the area occupied in 
the image by a given patch is affected by foreshortening. The actual surface area is 
proportional to 1 /(s» • v), where s,' is the normal of the patch, while v is the vector 
pointing towards the viewer (Figure 4). 

Measurements of surface orientation from images will not be perfect, since they 
are affected by the noise in brightness measurements. Similarly, surface orientations 
obtained from range data will be somewhat inaccurate. Consequently the impulses 
on the Gaussian sphere will be displaced a little from their true positions. The 
expected density on the Gaussian sphere will nevertheless tend to be equal to 
the inverse of the Gaussian curvature. One cannot, however, expect the impulses 
corresponding to a planar surface to be coincident. Instead, they will tend to form 
a small cluster. To be precise, the effect of noise is to smear out the information on 
the sphere. The extended Gaussian image is convolved with a smoothing function 
of width proportional to the magnitude of the noise. 
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Figure 13. Mapping of discrete patches on an object onto the Gaussian sphere. 
The patches in this case correspond to a regular tesselation of the image plane. 
Since the patches lie on a conical surface they contribute to the extended Gaussian 
image along a small circle. 


4.1. Using Object Models 

Extended Gaussian images also have to be computed for surfaces of prototypical 
object models. In this case it is best to find a convenient way to parameterize the 
surface and break it up into many small patches. Suppose the surface is given in 
terms of two parameters u and v as r(u, u). Then at the point (u,v) we see that r u 
and r„ are two tangents (Figure 14). The cross product of these tangents is normal 
to the surface. The unit normal 

r u X r v 

n = - 7 

|r u X it, I 

allows us to determine to which point on the Gaussian sphere this patch corresponds. 
Suppose that we divided the range of u into segments of size Su and the range of 
v into segments of size 6v. Then the area of the patch, 

6A — |r u X r„| Su6v, 

can be used to determine what contribution it makes to the corresponding place on 
the Gaussian sphere. Note that we do not have to explicitly compute the Gaussian 
curvature or take second partial derivatives. 
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Figure 14. A surface normal can be computed by taking the cross-product of 
two tangent vectors. The tangent vectors can be obtained by differentiation of the 
parametric form of the equation of the surface. 


5. Tesselation of the Gaussian sphere: Orientation Histograms 


It is useful to divide the sphere up into cells in order to represent the information 
on the Gaussian sphere in a computer. Ideally the cells should satisfy the following 
criteria: 

1. All cells should have the same area. 


2 . All cells should have the same shape. 

3. The cells should have regular shapes that are compact. 

4. The division should be fine enough to provide good angular resolution. 

5. For some rotations the cells should be brought into coincidence with themselves. 

Cells which are compact combine information only from surface patches which have 
nearly the same orientation. Elongated cells of the same area combine information 
from surface patches which have more widely differing orientations. The area of a 
regular polygon with n sides inscribed in a circle of radius r is 


7 rr 


2 


sin(27r/n) 

(27r/n) 


So the area of a hexagon incribed in a circle is (3\/3/2)r 2 , twice that of a triangle 
inscribed in the same circle. Tesselations with near-triangular cells will thus combine 
information from orientations which are \/2 times as far from the average as do 
tesselation using near-hexagonal cells. 
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Figure 15. The Gaussian sphere can be divided into cells along meridians and 
lines of longitude. The resulting cells do not have the same areas however and only 
align with each other for certain rotations about the axis through the poles. 


If cells occur in a regular pattern, the relationship of a cell to its neighbors will 
be the same for all cells. Such arrangements are to be preferred. Unfortunately, it 
is not possible then to simultaneously satisfy the criteria listed above. 

A simple tesselation consists of a division into latitude bands, each of which is 
then further divided along longitudinal strips (Figure 15). The cells could be made 
more nearly equal in area by having fewer at higher latitudes, or by making the 
latitude bands wider there, or both. One advantage of this scheme is that it makes 
it easy to compute to which cell a particular surface normal should be assigned. 
Still, this arrangement does not come close to satisfying the criteria stated above. 
In particular, the cells are brought into alignment only for a few rotations about 
the axis of the globe. Rotations about any other axis can not bring the cells into 
alignment. 

5.1. Tesselations Based on Regular Polyhedra 

Better tesselations may be found by projecting regular polyhedra onto the unit 
sphere after bringing their center to the center of the sphere [31] Regular polyhedra 
are uniform and have faces which are all of one kind of regular polygon (They are 
also called the Platonic solids). [8, 10, 19, 23, 25, 26, 30]. The vertices of a regular 
polyhedron are congruent. A division obtained by projecting a regular polyhedron 
has the desirable property that the resulting cells all have the same shape and area. 
Also, all cells have the same geometric relationship to their neighbors. In the case of 
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Figure 16. Tesselations of the Gaussian sphere using (a) the regular dodecahedron 
and (b) the regular icosahedron. 






the dodecahedron, the cells are even fairly well rounded. The dodecahedron, however, 
has only twelve cells (Figure 16a). Even the icosahedron, with twenty triangular 
cells, provides too coarse a sampling of orientations (Figure 16b). Furthermore, 
its cells are not well rounded. Unfortunately there are only five regular solids 
(tetrahedron, hexahedron, octahedron, dodecahedron, and icosahedron). 

One can go a little further by considering semi-regular polyhedra. A semi-regular 
polyhedron has regular polygons as faces, but the faces are not all of the same 
kind [8, 10, 19, 23, 25, 26, 30] (They are also called the Archimedean polyhedra). 
As for regular polyhedra, the vertices are congruent. There can be either two or 
three different types of faces and these have different areas. An illustration of a 
tesselation using a semi-regular polyhedron is provided by a soccer ball (Figure 
17a). It is based on the truncated icosahedron, a semi-regular polyhedron which 
has 12 pentagonal faces and 20 hexagonal faces. 

Unfortunately there are only 13 semi-regular polyhedra (The five truncated 
regular polyhedra, cuboctahedron, icosidodecahedron, snub cuboctahedron, snub 
icosidodecahedron, truncated cuboctahedron, rhombicuboctahedron, truncated 
icosidodecahedron, and the rhombicosidodecahedron). Overall, these objects do 
not provide us with fine enough tesselations. The snub icosadodecahedron has the 
largest number of faces, but each of its 80 triangles is much smaller than each of 
its 12 pentagons. 

The edges of a semi-regular polyhderon are all the same length. One consequence 
of this is that the different types of faces have different areas. The area of a regular 
polygon of n sides and edge-length e equals 


21 



Horn 


EGI 



Figure 17. Tesselation of the Gaussian sphere using (a) the truncated icosahedron 
and (b) the pentakis dodecahedron. 

ne 2 

4tan(7r/n)’ 

so it is very roughly proportional to n 2 . This is a problem generally with semi-regular 
polyhedra. It is sometimes possible to derive a new polyhedron which has the same 
adjacency relationships between faces as a given semi-regular polyhedron but also 
has faces of equal area. The shapes of some of these faces then are no longer regular, 
however. 

If we desire a finer subdivision still, we can consider splitting each face of 
a given tesselation further into triangular facets. If, for example, we split each 
pentagonal face of a dodecahedron into five equal triangles we obtain a pentakis 
dodecahedron with 60 faces (Figure 17b). This happens to be the dual of the 
truncated icosahedron, discussed above. If we apply this method instead to the 
truncated icosahedron we construct an object with 180 faces. This object, as well 
as the pentakis dodecahedron, form suitable bases for further subdivision, as we 
shall show later. 

To see how fine a division we might need, let us calculate the angular spread 
of surface normals which map into a particular cell. If there are n equal cells, then 
each one will have area 

A = (4?r )/n, 

since the total area of the unit sphere is 4ir (This area equals the solid angle of the 
cone formed by the cell when connected to the center of the sphere). The shape 
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which minimizes the angular spread for given surface area is the circular disc. The 
area of a circular disc on the unit sphere is 

A — 27r(l — cos#), 

where 6 is the half-angle of the cone formed by the disc when connected to the 
center of the sphere. If 9 is small, the area can be approximated by 

A !=» 7 r9 2 . 

Thus if there are many cells and if they could be made circular, the angular spread 
would be 

9 2 /\fn. 

The best we can hope for, however, are near-hexagonal cells. The area of a hexagon 
inscribed in a circle of radius r is (3\/3/2)r 2 , as already mentioned. The area of the 
circle, 7rr 2 , is about 20 % more. So a hexagonal shape has a spread which is 

= 1.0996... 

as large as that of a circular shape of equal area. A lower bound on the angular 
spread for a tesselation with n cells then is 



For n — 60, for example, the spread is at greater than 16.2 degrees. One should 
also remember that the spread for triangular cells is even more, namely y/2 times 
that for hexagonal cells. 

5.2. Geodesic Domes 

To proceed further, we can divide the triangular cells into four smaller triangles 
according to the well known geodesic dome constructions [19, 27, 31]. We attain high 
resolution by relenting on several of the criteria given above (Figure 18). Specifically, 
the cells of a geodesic tesselation do not all have the same area and shape. The cells 
are also not compact, being shaped like (irregular) triangles. The duals of geodesic 
domes are better in this respect, since they have facets that are mostly (irregular) 
hexagons, with a dozen (regular) pentagons thrown in. Tesselations of arbitrary 
fineness can be constructed in this fashion. The pentakis dodecahedron is a good 
starting point for a geodesic division, as is the object constructed earlier from the 
truncated icosahedron by dividing the faces into triangles. 

Each of the edges of the triangular cells of the original polyhedron are divided 
into / sections, where / is called the frequency of the geodesic division. The result 
is that each face is divided into f 2 (irregular) triangles. Tesselations where the 
frequency is a power of two are particularly well suited to the method suggested 
here, as we see next. 

One has to be able to efficiently compute to which cell a particular surface 
normal belongs. In the case of the tesselations derived from regular polyhedra, one 
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Figure 18. Tesselation of the Gaussian sphere using a frequency four geodesic 
tesselation based on the icosahedron (There are 16 X 20 = 320 faces). 


first computes the dot-product of the given unit vector and the vector to the center 
of each cell (These reference vectors correspond to the vertices of the dual of the 
original regular polyhedron). This gives one the cosine of the angle between the 
two. The closest reference vector is the one which gives the largest dot-product. 
The given vector is then assigned to the cell corresponding to that reference vector. 

In the case of a geodesic dome, it is possible to proceed hierarchically, 
particularly if the frequency is a power of two. The geodesic dome is based on some 
regular polyhedron. The appropriate facet of this polyhedron is found as above. 
Next, one determines into which of the triangles of the first division of this facet 
the given unit normal falls. This can be done by considering which dot-product has 
the second largest value. No new dot-products need to be computed. The process 
is then repeated with the four triangles into which this facet is divided, and so on. 
In practice, lookup table methods can be used, which, while not exact, are very 
quick. 

Let the area occupied by one of the cells on the Gaussian sphere be w (in 
the case of the icosahedron u — 47r/20). The expected number of surface normals 
mapped into a cell equals 

pu)\V\, 

for a convex object, where 73 is the average of G(£, rj) over the cell. 
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It is clear that the extended Gaussian image can be computed locally. One 
simply counts the number of surface normals that belong in each cell. The 
expression for the Gaussian curvature, on the other hand, includes first and second 
partial derivatives of the surface function. In practice, estimates of derivatives are 
unreliable, because of noise. It is important therefore that the extended Gaussian 
image can be computed without estimating the derivatives. 

The values in the cells can be thought of as an orientation histogram. It has 
recently be brought to my attention that this is analogous to a scheme used for 
histogramming directions of dendrites on neurons [5]. 

The result can be displayed graphically using normal vectors on each of the 
cells to represent the weight of the accumulated surface normals. A frequency 
two sub-division of the pentakis dodecahedron provides 240 cells, enough for most 
practical purposes (Figure 19). The angular spread in this case is about 11.5 degrees. 
An alternative way to present the extended Gaussian image graphically is by means 
of a grey-level image where brightness in each cell is proportional to the count. 
The surface of the Gaussian sphere may be projected stereographically instead of 
orthographically in order to preserve the shapes of the cells. Their areas will be 
distorted however. 

6. Solids of Revolution 

In the case of the surface of a solid of revolution, the Gaussian curvature is 
rather easy to determine. The solid of revolution can be produced by rotating a 
(planar) generating curve about an axis (Figure 20). Let the generating curve be 
specified by the perpendicular distance from the axis, r(s), given as a function of 
arc length s along the curve. Let 9 be the angle of rotation around the axis. Now 
consider the Gaussian sphere positioned so that its axis is aligned with the axis 
of the solid of revolution. Let £ be the longitude and rj be the latitude on the 
Gaussian sphere. 

We can let £ correspond to 6. That is, a point on the object produced when 
the generating curve has rotated through an angle 9 has a surface normal that lies 
on the Gaussian sphere at a point with longitude £ — B. 

6.1. Gaussian Curvature of Solid of Revolution 

Consider a small patch on the Gaussian sphere lying between £ and £-f-<5£ in 
longitude and between rj and r} + 6rj in latitude. Its area is 

cos rj6£6r}. 

We need only determine the area of the corresponding patch on the object. It is 

rSBSs, 

where 6s is the change in arc distance along the generating curve corresponding 
to the change Sr) in surface orientation. The Gaussian curvature is the limit of the 
ratio of the two areas as they tend to zero. That is, 
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Figure 19. Orientation Histogram collected on a geodesic dome derived from 
the pentakis dodecahedron (There are 12 X 5 X 4 = 240 faces). This is a discrete 
approximation of the extended Gaussian image. The length of the vector attached 
to the center of a cell is proportional to the number of surface normals on the 
surface of the original object which have orientations falling within the range of 
directions spanned by that cell. 


,, cosr)6£6ri cos n 6n cos r] dr] 

K = lira -- —— — = lim - 

«»-o r 66 6s 6y-*o r 6s r as 

since 6$ — 66. The curvature of the generating curve, kg, is just the rate of change 
of direction with arc length along it [11, 20, 24]. So 


and hence 



K __ k g cos rj 
r 

It is easy to see that (Figure 21), sinr? = — r s , where r s is the partial derivative of 
r with respect to s. Differentiating with respect to s we get 
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drj d . . 

cos, * = S (_rJ;= “ rs< ' 

and so we obtain the simple formula 



In the case of a sphere of radius R, for example, we have r = Rcos(s/R) for 
—(tt/ 2)R < s < -\-(n/2)R. Thus r ss — —(r/i? 2 ) and K = l/i? 2 . 

For some purposes it is more useful to express the radius r as a function of 
^ the distance along the axis, rather than as a function of arc length along the curve. 

Let the distance along the axis be denoted by z. It is easy to see that (Figure 21), 
tan 77 = — r z , and so, differentiating with respect to s, 


27 







Horn 


EGI 



Figure 21. This figure shows the relationships between infinitesimal increments in 
arc length along the curve, distance from the axis of rotation and distance along 
the axis of rotation. 


2 dr) / » 

sec 


dz 

~Tzz~J~) 

da 


where from the figure we see that cos i) = z s , so that 

kg cos rj = 

and finally 

K 


dr) 4 

— cos/? = —r zz cos r), 
as 


r ZZ 


since 


r(l + rj)= ’ 


sec 2 rj = 1 + r l 


6.2. Alternate Derivation of Gaussian Curvature of a Solid of Revolution (*) 

We first need to review Meusnier’s theorem [11, 20, 24]. Consider a normal 
section of a surface at a particular point. It is obtained by cutting the surface 
with one of the planes including the local normal (Figure 22). Suppose that the 
curvature of the curve in which the surface cuts this plane is k^. Now imagine 
tilting the plane away from the normal by an angle rj (using the local tangent as 
an axis to rotate about). The new plane will cut the surface in a curve with higher 
curvature. In fact, it can be shown that the new curve has curvature 
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Figure 22. The curvature of the curve obtained by cutting a surface using 
an inclined plane is greater than that obtained by cutting it using a plane which 
includes the surface normal. Meusnier’s theorem tells us that the ratio of the two 
curvatures is equal to the cosine of the angle between the two planes. 


kn / cos r/. 

It is easy to see this in the case of a sphere, since a plane including the center cuts 
the sphere in a great circle, while an inclined plane cuts it in a small circle of radius 
proportional to the cosine of the angle of inclination. 

Now, let us return to the surface of revolution. It is not hard to show that 
one of the principal curvatures at a point on the surface will correspond to a cut 
through the surface by a plane which includes the axis of revolution. The curve 
obtained in this way is just the generating curve of the solid of revolution. So one 
of the two principal curvatures is equal to the curvature kq of the generating curve 
at the corresponding point. 


i 
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Now consider a plane perpendicular to the axis of revolution through the 
same surface point (Figure 23), It cuts the surface in a circle. The curvature in 
this plane equals (1/r), where r is the radius of the solid of revolution at that 
point. This horizontal plane, however, is not a normal section. Suppose that the 
normal makes an angle rj relative to this plane (The local tangent plane also makes 
an angle rj relative to the axis of revolution). Now construct the plane including 
the local normal which intersects the horizontal plane in a line perpendicular to 
the axis. This plane will be inclined r] relative to the one we have just studied. 
It also produces the second principal normal section sought after. By Meusnier’s 
theorem we see that the curvature of the curve found in this normal section is 
Kfir = (1/r) cosr]. Finally, the Gaussian curvature is found by multiplication to be 

K _ K G cos V 
r 

In the case of a sphere of radius R, for example, we have r — R cos rj and kq — 1/ R, 
so that K — 1/R 2 , as expected. 

To make this result more usable, erect a coordinate system with the 2 -axis 
aligned with the axis of revolution. The generating curve is given as r{z). Let 
the first and second derivatives of r with respect to z be denoted by r z and r zz 
respectively. It is easy to see that (Figure 21), tan 7] — r z , so that 

1 

COST] — 

Furthermore, 

__ r zz 

K ° ~ _ (1 +rj)3/*’ 

so that finally, 

jy- _ 

r(l + r2) 2 ‘ 

In order to use this result in deriving extended Gaussian images it is necessary to 
identify points on the surface with points on the Gaussian sphere. Suppose that we 
introduce a polar angle 9 such that e 

x — r cos 9 and y = r sin 9. 

Then a unit normal to the surface is given by 

(cos 9, sin 9, — r z ) T 

Equating this to the unit normal on the Gaussian sphere, 

(cos £ cos r], sin £ cos rj, sin rj) T , 

we get 




£ = 9 and tan?; = — r z . 
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Figure 23. If the solid of revolution is cut by a plane perpendicular to the axis of 
rotation, a circle is obtained. The curvature in this plane is just the inverse of the 
distance of the surface from the axis. The curvature of the corresponding normal 
section can be obtained using Meusnier’s theorem. 



6.3. Extended Gaussian Image of a Torus 

As an illustration we will now determine the extended Gaussian image of a § 
torus. Let the torus have major axis R and minor axis p (Figure 24). A point on 
the surface can be identified by 6 and s, where 0 is the angle around the axis of 
the torus, while s is the arc length along the surface measured from the plane of 
symmetry. Then, 

r = R -f- pcos(s/p), 

and 

r ss = —(1 /p) cos {s/p), 

so that 
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Figure 24. A torus obtained by spinning a circle around an axis. The resulting 
object is not convex. Its extended Gaussian image can be computed nevertheless. 


K _ r ss __ 1 cos {s/p) 

r p Rp cos(s/p)' 

Two points, P and P 1 (Figure 24), separated by 7r in 9, have the same surface 
orientation on the torus. The surface normal at one of these places points away 
from the axis of rotation, while it points towards the axis at the other place. 
Accordingly, two points on the object, 


(0»*) = (£,P»7) and {9,s) = (£ + n, p(n — ri)), 


correspond to the point (£, rj ) on the Gaussian sphere. The curvatures at these two 
points have opposite signs 


K = | 1 C0ST ? 

p R -f p cos rj 


and 


1 cost; 
p R — p cos 77 


The torus is not a convex object, so more than one point on its surface contributes 
to a given point on the extended Gaussian image. If we add the absolute values of 
the inverses of the curvature we get 


G(t7,£) = 


1 1 
K + Kl 


= 2Rp sec rj. 
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If we had added the inverses algebraically instead we would have obtained 


1 




which is twice the result for a sphere of radius p. 


The same results could have been found using 

«GC0s?7 

k = ~t~- 

since k g — —1 /p and r = R i pcos r), so that 


_ i cos rj n 

P [R ± P cos ri)' 

The extended Gaussian image of a torus has singularities at the poles. These 
correspond to the two rings on which the torus would rest if it were dropped onto 
a plane. All of the points on one of these rings have the same surface orientation. 


We can think of the Gaussian sphere as covered by two sheets of a Riemann 
surface, one corresponding to the inner half of the torus, closer to its axis of 
symmetry, the other corresponding to the outer half. The two sheets are connected 
to one another at the poles, branch points corresponding to the two rings mentioned 
above. There the Gaussian curvature changes sign. 

We may also note at this pdint that all tori with the same surface area, [Air^pR), 
have the same extended Gaussian image. 


6.4. The Unique Convex Object with (7(£, rj) = 2 sec 77 (*) 


While all tori with surface area 47 r 2 have the same extended Gaussian image, 

G[£,v) = 2 sec 77 , 

there is only one convex object which has that extended Gaussian image. It is a 
solid of revolution since <_?(£, 77 ) is independent of £. So we have on the one hand 

K = 1/2 cos 77 , 

and on the other hand 


so that 


K = *G cos 77 
r 


k g = r/2. 

The equation states that the curvature of the generating curves varies linearly with 
the distance from the axis of rotation. This deceptively simple equation represents 
a non-linear second order differential equation for r in terms of z since 

K _ r ** 

(l + r*)»/»' 

so that 
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and 


so that 


' zz 




d 1 __ r z r zz 

dz u + r *) 3/2 


d r 2 rr z 

dz A 2 ’ 


d 1 _ dr 2 

S v/TTi =S7 


or 

1 r 2 -f- c 2 

V 1 + r l 4 

where c 2 is a constant of integration. We now have reduced the problem to a 
non-linear first order differential equation for r in terms of z. If the object is to be 
convex and smooth at its poles, we expect r z —► oo as r —► 0. Thus c = 0. Next 
note that the term on the left equals cos r). So we also have 

r 2 

coarj — — 

4 

or, using an earlier expression for kq, 


kq = —a/cos rj. 

This is an implicit equation for the curve of least energy [14]! The curve of least 
energy is the curve which minimizes the integral of the square of the curvature Ke¬ 
lt can be solved for z in terms of r to yield 

z = V2[2E(cos~ 1 (r/2), 1 / \/2) — F(cos —1 (r/2), 1/V^)]. 

where E and F are incomplete elliptic integrals. If we let s be the arc length along 
the curve we can also write the solution in Whewell form 

s — V2F(cos~~ 1 \J cos Tj, l/\/2), 


or Cesaro form 

s = \/2.F(cos~ x (— kq), l/s/2). 
The length of the curve from the pole to the equator is 


\f2K{l/\/2) — >/2.ftT(sin(7r/4)) — 


r(i/4) 2 

2\/27r 


where K is the complete elliptic integral of the first kind, and T is the gamma- 
function [20]. The height from equator to the pole is 
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Figure 25. The unique convex object with the same extended Gaussian image 
as a torus has an interesting shape. It is a solid of revolution whose generating curve 
is the curve of least energy. This is the shape which a uniform bar constrained to 
pass through two points in space with given orientation will adopt. 



W = 


(2tt) 3 / 2 

F(l/4) 2 


while the maximum radius is 

H = 2. 

The minimum radius of curvature equals one, so that a circle tangent at the 
outermost point is also tangent at the origin [14]. This circle, when rotated about 
the vertical axis, produces a torus with the same extended Gaussian image (Figure 
25). Both objects have total surface area 47T 2 . 


7. Gaussian Curvature in the General Case 

When the object is not a solid of revolution we need to work a little harder 
to obtain the Gaussian curvature. Let x — x(u, v ), y — y(u, v), and z — z(u, v) be 
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parametric equations for points on a given surface. Let r = (x,y,z) T be a vector 
to a point on the surface. Then 

dr ' dr 

r« = — and r„ = — 

du dv 

are two tangents to the surface, as already noted earlier. The cross product of these 
two vectors, 

n = r u X r w> 

will be perpendicular to the local tangent plane (Figure 14). The length of this 
normal vector squared equals 

n 2 = n • n = (r u • r u )(r v • r„) — (r u • r„) 2 

since (a X b) • (c X d) = (a • c)(b • d) — (a • d)(b • c). A unit vector n = n/n can be 
computed using this result. 


7.1. Gaussian Curvature from Variation in Normals (*) 

The Gaussian curvature is the limit of the ratio of the area of a patch on the 
Gaussian sphere to the area of the corresponding patch on the surface, as the area 
shrinks to zero. Consider an infinitesimal triangle formed by the three points on 
the surface corresponding to (u, v ), (u -j- 6u,v), and ( u , v -J- <$t>). The lengths of two 
sides of this triangle are 

jr u |<5u and |r„|<5v 

while the sine of the angle between these sides equals 

K X r w | 

M W 

so that an outward normal with size equal to the area of the triangle is given by 

~(r u X r v )6u6v — -nSuSv. 

2 ' 2 

To determine the area of the corresponding triangular patch on the Gaussian sphere 
we need to find the unit surface normals at the three points. The unit surface 
normals will be 

n, and 

if we ignore terms of higher order in Su and 6v. Here n u and are the partial 
derivatives of n with respect to u and v. Note that h u and n v are perpendicular to 
n. The area of the patch on the Gaussian sphere equals the magnitude of 

1 

-(n u X n„)<5u6v 
z 

by reasoning similar to that used in determining the area of the original patch on 
the given surface. We need to find n u and n v to compute this area. Now 

d n nn u — n n u 
dun n l 
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I 




From n 2 = n • n we get 
so that 

fi u = 
and 

n« = 


nn u — n • n u , 

(n • n)n u — (n • n u )n (n X n u ) X n 






(n • n)n„ — (n • n„)n (n X n v ) X n 


n° 


n J 


since (a X b) X c = (a • c)b — (b • c)a. Then 


n 


X n„ = -g [(n • n)(n u X n„) + (n • n u )(n v X n) + (n * n„)(n X n u )] 


n 


or 


n n Xn v — -4-[nn u n v ]n 
n q 

since [a b c]p — (a • p)(b X c) -j- (b • p)(c X a) + (c • p)(a X b), where [a b c] = (a X b) • c. 

This shows that the patch on the Gaussian sphere has the same orientation as 
the patch on the surface, as it should. An outward pointing normal of size equal to 
the area is given by 

1 1 


2 n 4 


n n u n„]n<5u6u. 


The ratio of the two areas, the Gaussian curvature, is then just 

[n n« n„] 


K 


n ’ 


Now 
so that 


n =.r tt X r«, 

n u = r uu X r„ + r u X r uv , 


and 


Hti — Tuv X -f- r„ X Tvv 

Now using (aXb)x(cXd) = [abd]c—[abc]d or(aXb)x(cXd) = [acdjb—[bcd]a 
we get 

DuX^t/ — —[ r uu r v r vl i]r^ -j- [r uu r w ]r u [r uu r^ Ptx.JPw ~ l - [^u iw T* [r« ^uv rvy]r u 

so that 

[n n u n v ] = n • (n u X n^) === [rur^ruuIJryrvrt,^] [rur v r U t,] , 

and finally 


tj- _ [ r u iy r MU ] [r u t v r w ] [ru r v r ut ,] 2 

|r« X r v | 4 

This result can be used to derive the expression for curvature of a solid of revolution 
in a more rigorous fashion. 
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7.2. Fundamental Forms of a Surface (*) 


Let, as before, 

»_ fu X r® 

|r u X r w | 

be the unit surface normal vector. The first fundamental form of a surface gives 
the square of the element of distance as [20] 

ds 2 = |dr| 2 = E(u, v ) du 2 -f 2 F(u, u) du dv -f- G(u, v) dv 2 . 

The second fundamental form of a surface gives the normal curvature using the 
equation [20] 

—dr • dn — L(u, v ) du 2 -f- 2 M(u, u) du dv -f- N(u, v) dv 2 . 

The coefficents can be expressed in terms of derivatives of r as follows, 

E — r u - r u , F — t u - r v , and G = r v -t v 


and 


or, 


L = r u -n u , M ~ r u • n v = t v ■ n u , and N = r v -n v 


M _ Jr^jyr^L and jV _. Kr,r„j 
, ' / — —■ — ' 


\/EG — F 2 


so that [20] 


K 


Veg — f^ 
LN — M 2 


VEG - F 2 


EG — F 2 * 

Finally, if the surface is given as z(x,y), the above reduces to the familiar 

»2 


K 


z xx z yy z xy 


(l — 1“ Z \ ”4“ Z y^j 


2N2 ’ 


7.3. Application of the General Formula to the Ellipsoid (*) 

In the case of the ellipsoid we have, as discussed before, 

r = (a cos 6 cos <j>, b sin 6 cos (f> , c sin <j >) r , 
to = (—a sin 6 cos <f), b cos 0 cos <f>, 0) T , 
t^ — (—a cos 6 sin </>, —6 sin 9 sin <fi, c cos <j>) T , 
and 

roe — (—cl cos 9 cos (j), —b sin 9 cos <j>, 0) T , 

t$ ip = (a sin 9 sin <f>, —b cos 9 sin <f>, Oj^, 

rpp — (—a cos 9 cos (j), —b sin 9 cos <j>, —c sin <j>) T . 

A surface normal can be found by taking cross products, 

n = to X Tp — (be cos 9 cos <f>, ca sin 9 cos <j>, ab sin <f>) T cos <f >, 
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and the coefficients of the first fundamental form are, 


E = tq ■ r$ = (a 2 sin 2 9 -f- b 2 cos 2 0)cos 2 <j >, 

F = tq • — (a 2 -f- 6 2 ) sin 9 cos 9 sin <f> cos <f>, 

G == = (a 2 cos 2 6 + b 2 sin 2 9) sin 2 4> -f- c 2 cos 2 <ft. 

Hence 

EG — F 2 — [(be cos 9 cos <ft) 2 -f- (ca sin 9 cos (ft) 2 -f- ( ab sin <f >) 2 ] cos 2 <j>. 
To compute the coefficients of the second fundamental form we need, 

[re r^ ree] = n • rg$ = —abc cos 3 <j>, 

[r^ rgft\ = n • r^ — 0, 

[re r^, r^^j — n * r^^ — abc cos 


and so 


Finally then, 


[re i> ree][re r^, r^] — [re r^ r^] 2 = (a6c cos 2 (ft) 2 . 


LN — M 2 = 


(abc cos 2 <ft) 2 

EG-^W' 


and 


LN — M 2 
EG — F 2 


abc 

(be cos 9 cos (ft) 2 + (ca sin 9 cos (ft) 2 -f- (ab sin <ft) 2 


This result was used earlier in the discussion of the extended Gaussian image of 
the ellipsoid. 
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